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INTRODUCTION 


Our ability to predict the performance of Army systems in natural 
and disturbed environments depends on our understanding of the chemistry 
and physics of the atmosphere. In order to improve our knowledge in this 
area, the BRL undertook the Combined Stratospheric Measuring Program 
(COSMEP) in the summer of 1972. The goal of this program is to make 
near simultaneous measurements of stratospheric minor constituents, both 
charged and neutral. This report deals with one of the techniques used 
to obtain mixing ratio profiles for minor neutral constituents. 


The use of infrared remote sensing techniques for studying atmos- 
pheric composition has increased greatly within the last decade. One of 
these techniques is to analyze transmission spectra of the earth's atmos- 
phere which are made using the setting sun as a source. This technique 
has the advantages that 1) the long optical paths through the atmosphere 
increase the sensitivity of this technique to minor constituents, 2) the 
long paths through the atmosphere decrease the sensitivity of this 
technique to localized contamination, especially water vapor, and 3) the 
geometry associated with this technique allows the vertical mixing ratio 
profiles for minor constituents to be established.’ 


In order to use this technique to derive mixing ratio profiles, one 
must know the quantity of gas in the optical path, the so-called atr 
mass. Since this quantity 1s so fundamental, it has been calculated by 
many investigators for a variety of conditions, under a variety of dif- 
ferent assumptions. 


One assumption which greatly simplifies the calculation of the air 
mass is that atmgspheric refraction is negligible. In this case, the 
Chapman function can be used to calculate the air mass, and simple 
geometry can be used to compute the minimum altitude for a grazing ray. 


The subject of atmospheric refraction has been treated by many 
authors and was originally discussed by Laplace. It has been reviewed 
recently by Kondratyev, and will not be reviewed here, except to the 
extent that it contributes to the air mass calculations. 


' Drayson, S. R., Bartman, F. L., Kuhn, W. R. and Tallamraju, R., 
"Feasibility of Satellite Measurements of Stratospheric Minor Con- 
Stituents by Solar Occultation,"' Technical Report 011023-2-T, High 
Altitude Engineering Laboratory, University of Michigan, Oct. 1975. 


*Farmer, C. B., "Infrared Measurements of Stratospheric Composition," 
Can. J. Chem., 32, 1544 (1974). 


chapman, S., 'The Absorption and Dissociative or Ionizing Effect of 
Monochromatic Radiation in an Atmosphere on a Rotating Earth. MII. 


Grazing Incidence,"Proc. Phys. Soc. (London) 43, 483 (1931). 


“vondeatyey, K. Ya., "Radiation in the Atmosphere,'', Academic Press, 
1969. 


One of the earliest tabulations of,air mass which took refraction 
into account was performed by Bemporad,~’ and has been generally accepted 
by astronomers. This calculation was performed for a ground based 
observer and for zenith angles less than or equal to 90°. A closed form 
expression, approximating Bemporad's results to within 1.5% for zenith 
angles less than or equal to 87.5°, has been given by Nagel. 


Since air mass and refraction depend on the meteorological condi- 
tions of pressure, temperature, and humidity, as well as on wavelength, 
one would expect that the computed air mass and refra tion would vary 
depending on the assumed atmospheric profile. Kasten ’° has yyudied the 
etfect of density variations on relative air masses. Link ’  ’ ~ 
has published a large number of computations, the most comprehensive of 
which ~ tabulates the refraction, air mass, and related quantities for a 
variety of reference atmospheres, observer and source heights, and 
zenith angles (70 to 90°). The index of refraction used in the above 
calculations corresponds to light with a wavelength of 500 nm. Formulae 
for correcting the tabulated refraction are given for wavelengths 
between 400 and 700 nm. 


None of the tables described above gives the refraction and the air 
mass for an observer above the surface and for zenith angles greater 
than 90°, the so-called oceceultation geometry. Several investigators 
have filled this gap by developing computer programs to calculate refrac- 
tion and air mass under specific conditions. 


Cc 
~Bemporad, A., "Search for a New Empirical Formula for the Represen- 


tation of the Variation of the [ntensity of Solar Radiation with 
Zenith Angle;"Rend. Acc. Lincei, 6, 66 (1907). 


Ge. ; : ; F sae le 
List, R. J.,"Smithsonian Meteorological Tables,"' 6th Revised Edition, 


1949, Fifth Reprint, 1971, Smithsonian Publication 4014. 


‘Nagle, M. R., ‘Improved Approximation of Bemporad's Atr Mass Func- 


tion," Appl. Opt. 13, 1009 (1974). 


Kasten, F., 'A New Table and Approximation Formula for the Relative 


Optical Airmass,'' Arch. Meteorol. Geophys. u. Bioklimatol. B14, 
206 (1966). 


9 ae ; ; 
“Kasten, F., "Effect of Variation of the Vertical Air Density Profile 


on the Relative Optical Air Mass," Arch. Meteorol. Geophys. u. 
Bioklimatol. B15, 62 (1967). 


{) 
mitt. D F. and Sekera, Z., "Dioptrische Tafeln der Erdatmosphare," 
Publikace Prazske, Hvezdarny Prague (1940) 


wo F., "Extension des Tables Dioptriques De L'Atmosphere 


Terrestre," Prometheus, Prague, 1947. 


9 
kt os F. and Neuzil, L., "Tables of Light Trajectories in the Terres- 


trial Atmosphere," Hermann, Paris, 1969. 
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Selby and McClatchey** have included refraction in their multi- 
purpose computer code LOWTRAN 2. This code is flexible enough to include 
almost any remote sensing geometry. They have defined six reference atmos- 
pheres with respect to pressure, temperature, and water vapor content. 
These reference atmospheres are used to compute the index of refraction 
as a function of altitude. They also took into account the deviation of 
the earth's surface from that of a sphere. 


Watson and Yin?4 have studied the effects of refraction on solar 
occultation measurement s of atmospheric gases. They have included anom- 
alous dispersion, and have computed the "infinite resolution" and instru- 
mental line shapes for three H,O lines under a variety of atmospheric 
models. They formed the absorption coefficient in a line-by-line and 
layer-by-layer process, and therefore did not use a single layer approx- 
imation of the Curtis-Godson type. 

16,17 18 ' 

Hymus and Lloyd, and Miller have computed the refraction and 
extinction of a near ultraviolet tangential ray. Since different 
processes govern the extinction of infrared and ultraviolet radiation, 
their gesylts are not directly applicable to the IR remote sensing problem. 
Treve ~* has described a ray-tracing technique which can be used to 
compute the atmospheric refraction and air mass for a number of 

different source-observer geometries, and for a wavelength range of 0.2 

to 20 micrometers. 


1S cntie, J. E. A., and McClatchey, R. A., "Atmospheric Transmittance 


from 0.25 to 28.5 um: Computer Code LOWTRAN 2,'' Technical Report 
AFCRL-72-0745, Environmental Research Paper No. 427, 29 Dec 1972. 


'4watson, J. K. G. and Yin, P. K. L., "The Effects of Refraction and Dis- 
persion on High Altitude Measurements of Atmospheric Gases," 
Interim Report R.F. Project 3460, The Ohio State University Research 
Foundation, October 1973. 


> 

Calfee, Robert F., "Anomalous Dispersion Calculated for Atmospheric 
Water Vapor," Appl. Opt. 7, 1652 (1968). 

16 


Hymus, F. C., and Lloyd, K. H., "Extinction of a Tangential Ray, by 
Ray Tracing Through The Atmosphere,'' Technical Note HSA-152, Aus- 
tralian Defence Scientific Service, Weapons Research Establishment, 
May 1969. 

eine K. H., "Absorption of a grazing Ray Calculated by Ray Trac- 


ing,'' Planetary and Space Sciences 17, 1683 (1969). 


1S tee ND. E., "Stratospheric Attenuation in the Near Ultraviolet," 


Proc, Roy. Soc. A 301, 57 (1967). 


Ie ere, Y. M., "New Values of the Optical Air Mass and the Refraction 
and Comparison with Previous Tables,"' Technical Report ESD-TDR-64-105, 


Mitre Corp. 1964. 

2 

Ow ies Y. M., "A New Method for the Calculation of the Atmospheric 
Refraction and the Optical Air Mass," Unpublished. 
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This report describes a series of air mass and refraction calcula- 
tions, performed using a modified and extended version of Treve's ray- 
tracing technique. In an attempt to assess the accuracy of this technique 
several checks and comparisons were performed, and are also described in 
this report. Finally, tables giving the refraction, air mass, and 
Curtis-Godson effective temperature and pressure for both a summer and 
winter reference atmosphere are given in Appendix A. Tables in Appendix 
A cover an observer altitude range of 10 to 50 km in steps of 1 km. For 
each altitude, the zenith angles go from 80° to 89° in steps of 1.0°, 
and from 89° to the disc in steps of 0.2°. A statement listing for the 
computer program used to generate the tables in Appendix A is given in 
Appendix B. 


THEORY 


The technique used for these calculations is basically that of 
Treve, with several modifications. A brief description of that technique 
will be given here with a somewhat more detailed description of the 
modifications, , Bar more details, the reader is referred to the original 
descriptions. ~’” 


The model assumes a dry atmosphere, which is made up of spherically 
symmetric concentric shells, of known temperature and pressure. The 
atmosphere is assumed to be an ideal gas* in hydrostatic equilibrium. 
Together with the above assumptions, Snell's law can be used to derive 
the relation (see reference 13): 


Nn, Fr. $in Zz. = constant Ci) 
., Y 1 
where: n, 1s the phase refractive index of level i. 
7 1s the distance from the center of the earth to level 1. 
and a is the angle of incidence with respect to the local 
normal. 


The phase refractive index for dry air is given by reference 22: 
(n, - 1.0) x 10° = 64.328 + 29498.1 (146 - acy? 

2,-1 (2) 

+ 255.4 (41 - 1/)d°) 5 

for a temperature of 288.16K and a pressure of 101325.0 N/m”. The 

wavelength A is given in micrometers. This expression for the index of 

refraction neglects the effects of anomalous dispersion as well as the 

contribution of water vapor. Anoma lous dispersion, has been studied by 

Calfee, and Watson and Yin, and has been shown to be negligible 

for stratospheric conditions. The corre¢tion to the index of refractions 

for water vapor has been given by Owens and has also been determined 

to be small for stratospheric conditions. 


Strictly speaking, this is not true (see Ref. 21), but to the degree of 
accuracy that we are interested in here, this approximation is a good one. 


Owens, James C., "Optical Refractive Index of Air: Dependence on 
Pressure, Temperature, and Composition," Appl. Opt. 6, 51 (1967). 


ae 
“" Edlen, B., "The Index of Refraction for Dry Air in the Visible and 
Infrared,"' Metrologia 2, 71 (1966). 
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The index of refraction for each shell was computed from Equation (2) 
assuming it to be a function of density alone. This dependence is given 
by the Lorentz-Lorenz relation: 


ee ae 
— (8+) = constant (3) 
n + 2 


The density p in this expression is obtained from the perfect gas law; 


M P 
ye ET re 


where: M is the mean molecular weight of air = 28.966, 
R* is the universal gas constant = 8314.39 J/K/kg-mole, 
P is the pressure in Ty ; 
and T is the temperature in K. 
For the calculations performed in this study, various temperature 
profiles were used and the corresponding pressure profiles were computed 


from the following relations derived [assuming an ideal gas - Reference (23) ] 
from the hydrostatic equation: 


N 
‘ Fon 
( O ) R*L 
= -. + L(H - H,) L # 0 
» (H - H)) 
aa gil Gel: a aa a 
where: P is the pressure at the base level, 


H is the geopotential altitude, 
H is the geopotential altitude of base, 
L is the gradient of temperature T with geopotential altitude H, 


and g is the acceleration of gravity - assumed to be a function 
of latitude, but not of altitude. 


Since the conversion factors between geopotential and geometric alti- 
tude as well as the acceleration of gravity g, are functions of latitude, 
Treve's model was extended to, cgyer all latitudes by using the appropriate 
expressions given by Lambert: ’ 


+7 6, Standard Atmospheric Supplements, 1966. 
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rH (6) 


4. {Hh} =f 
gM) = Ty8, ‘. 
G 
Yr. « g 
+ = 4 ° _9 

H, (2) se ane f (7) 

> 
where: A, = altitude in geopotential meters (m'), 
an = altitude in geometric meters (m), 


2 “ - 
G = 9.80665(m)*~ (sec) hw) ‘ 


-2 
= sea level acceleration of gravity at latitude $, (m)(sec) ~, 


Bo 
Ty = effective earth radius for latitude 6, (m), 
g, = 9.806160 {1.0 - 0.0026373 cos 2¢ 
+ 0.0000059 con® 2¢ | (8) 
and 2 Be 
r, = Oe ee 9) 


6 


% 3.085462 x 10°°+ 2.27 x 107? GOS 26-- < xX 1° +4 cos 4¢ 


This expression for the effective earth radius r, is not in general equal 
to the radius of the earth at the given latitude: It is rather a ficti- 
cious quantity, designed to satisfy boundary conditions and to take into 
account the non-spherical aie 3 of the earth, its density distribution, 
and centripetal acceleration. 4 


The geocentric radius vector to the surface of the ellipsoid, which 


approximates the surface of the earth, is computed from the relationship: “4 
Te Ge ee 
a cos > + Sin 
R($} . = “—— pea (10) 
a’ cos » + b” sin’ ¢ 


where: a equatorial radius = 6378.160 km 


b = polar radius 6356.775 km 


and t latitude 


fi 
24 ankey, John R., "The American Ephemeris and Nautical Almanac 1974," 


U.S. Government Printing Office, 1971. 
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The earth is here approximated as a sphere, with an effective radius 
R(¢), which depends on latitude. Table 1 shows a comparison of T 5 and 
R(¢) as a function of latitude. 


TABLE 1. COMPARISON OF ry AND R(¢) 


Latitude R(d) km ty km 
0.0 6378.160 6334.984 
30.0 6372. 847 6337.838 
45.0 6367.512 6356. 360 
60.0 6362.155 6367.103 
2.0 6356.775 6377.862 


For the purposes of ray tracing, it 1s important to use the geometric 
radius vector R(¢) rather than the effective radius ry. For altitudes 
of less than 100 km, the error in geopotential altitude H will be less 
than 0.01 km. 


The relevant geometry for this type of calculation is shown in 
Figures 1 and 2. H., is the altitude of the observer in geometric meters 
above the earth's surface. The astronomical zenith angle (z) is the angle 
between the local normal and an unrefracted ray from the source. The 
apparent zenith angle (z') is the angle of arrival of a refracted ray 
from the source. The refraction is the difference (z - z'). 


For the case that z' is greater than 90° (Figure 2) the grazing 
altitude (minimum height) is less than the observer height. As long as 
the density profile is monotonic decreasing with altitude, the minimum 
height for a refracted ray will always be greater than that for an unre- 
fracted ray. From Figure 2 and simple geometric consideration, the 
minimum altitude for an unrefracted ray can be seen to be: 


Han = (R(o) + Hop? sin (mt - z) - R(¢) ere 


Thus the unrefracted minimum height is a function of observer height, 
zenith angle, and earth radius (which is a function of latitude). 


The air mass, or the quantity of atmosphere traversed by the ray, 
is defined as: 


Source 
air mass = aes a 29 N ds (12) 
24153550 % 10 
Observer 


where ds is the arcjength along the refracted ray, and N is the number 
density (molecule/m~),. Since the source is outside the atmosphere, inte- 
gration is cut off at the height of the highest level of the assumed 
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Sun 


Figure |. Refracted and unrefracted rays for zenith angles less than 99°. 
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Sun 


Figure 2. Refracted and unrefracted rays for zenith angles areater 
than 90°. 


protile, For pygposes of comparison, one standard air mass is here defined* 
as 2.15335 x 10°° molecules/m’. This corresponds to the numerical inte- 
gration of Equation (12) with z = 0 between 0 and 100 km for the 1962 U.S. 
Standard Atmosphere. Thus the present calculation gives an absolute air 
mass, to which other profiles can be compared, rather than a relative 

air mass, as originally computed by Treve. 


For the case when refraction is neglected, the air mass can be com- 
puted for large zenith angles from the Chapman function.~ Simplified,, 
approximations to this function have been derived by Smith and Smith, 
from expressions given by Fitzmaurice and modified for a non-exponential 
atmosphere. They are given by: 


air mass = Ue oy SS Ch(X_,x.) = fi2) 
2.15335 x 10°" pF 


n_ H 2 
59 [ (1/2) x! oY  eebety) Xp < 90° 


2,13533c5 % 10° 


5 H 4 
air mass = [(1/2)X ji/e occreeaeingaes: [2 n -n git ) -epfety)] (14) 
. 2.15335 x 10 
; ° 
ky + SY 
Where a, = number density at observer height, 
ny = number density at unrefracted grazing height = fee. 253 
H = scale height at observer altitude, 
i = scale height at grazing altitude, 
‘. = (R(o) + Fe) 2 
‘, = (R(o) + Reet! 
Xp = astronomical zenith angle at observer, 
Lf2 
= (X_/2) cos 
Rah, Jcos x, | 
and erfc(y) = 1.0 - erf(y) 


ee 


Note: The term air mass does not have a standard definition. Before 
comparing the values in this report to other work, the definitions must 
be compatible. 


oP osith. F. L. III, and Smith, Cody, "Numerical Evaluation of Chapman's 


Grazing Incidence Integral ch(X,x),'"" J. Geophys. Res. 77, 3592 (1972). 


eel eamaurice, John A., "Simplification of the Chapman Function for 


Atmospheric Attenuation," App. Opt. 3, 640 (1964). 
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These approximations take into account the fact that the air mass is most 
sensitive to the scale height and.gensity near the observer for zenith 


angles less than or equal to 90°. For zenith angles greater than 90°, 
the az mass 1s most sensitive to the scale height near the grazing alti- 
tude. For X = 800 this approximation is expected to be better than 1%. 


In order to compute the infrared absorption for a non-homogeneous 
ray path, such as used in this study, one must evaluate the absorption 
coefficient at each wavelength for each layer. This turns out to be a 
very time consuming computation. 


An alternative to this procedure is to define an ''equivalent" homo- 
geneous path at an effective pressure P(eff) and effective temperature 
T(eff), such that the calculated absorption is the same as that for the 
non-homogeneous case. For a uniformly mixed gas such as CO,, Sasamori 
has suggested defining these quantities as: 


Source / Source 
P dm dm 


P(eff) = 
Observer | Observer 
(13) 
Source / Source 
T(eff) = T dm / dm 


Observer / Observer 


Source | 
Where if dm gives the total integrated absorber mass, along the 
ray path. Thus the effective temperature* and pressure are essentially 
mass weighted average values. 


For the case that the gas in question is not uniformly distributed, 
the varying mixing ratios should be included inside the integrals of 
Equation (15). However, for the case that Z > 90°, the air mass increases 
rapidly with a small change in zenith angle. Thus for this case, the 
effective pressure and temperature are most sensitive to the mixing ratio 
at the tangent height, since most of the air mass is located within a 
few kilometers of the tangent height. 


* 
Strictly speaking, the,gffective temperature depends on the intensities 


of the spectral lines. For stratospheric occultation geometries, 

the temperature variation along a tangent ray will not be large. In 
this case, the mass weighted average temperature will be a good approx- 
imation to the effective temperature. 


“) swider, William Jr., ''The Determination of the Optical Depth at 


Large Solar Zenith Distances", Planet. Space Sci. 12, 761 (1964). 


28 osamori, T., "The Temperature Effect on the Absorption of the 15u 


CO, Band,” Sci. Rept. Tohoku Univ. Fifth Series 11, No. 3. 


mast t G. N., "Spectral Band Absorptance for Atmospheric Slant Paths," 


Appl. Opt. 2, 515 (1963). 
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VERIFICATION OF MODEL 


In order to verify the numerical integration of the ray tracing tech- 
nique and the above approximations to the Chapman function, the following 
calculation was performed. An isothermal exponential atmosphere 
(T = 272.344°K) was chosen so that X_ would be 800.00 for an observer 
at 10 km. A ray trace was performed?on this atmosphere ignoring the 
effects of refraction (n. = 1.000000000000). The Chapman function derived 
from the numerically integrated ray traced air mass is compared to the 
value calculated by Wilkes.” The value for the Smith and Smith approxi- 
mation is also compared to the Wilkes value. Table 2 shows the results 
of these comparisons. It can be seen that the Smith and Smith approxi- 
mation agrees well with the Wilkes value, which is good to three decimal] 
places. In all cases, the agreement is better than 1/2%. The ray 
tracing technique gives better results for larger zenith angles, but in 
all cases, the difference is less than 1/2%. This agreement gives an 
estimate of the error associated with the numerical techniques used to 
approximate the ray trace air mass integral. 


In an attempt to verify the accuracy of this technique for a non- 
exponential atmosphere, the air mass for a ground based observer was com- 
puted for a series of zenith angles, assuming the 1962 U.S. Standard 
Atmosphere (58 layers between 0 and 100 km). For comparison, the same 
calculation was performed for the case of no refraction. The results 
are tabulated in Table 3. There it can be seen that inclusion of refrac- 
tive effects can increase the computed air mass by up to about 10% for 
z = 90°. This is probably an upper limit since water vapor, which was 
neglected in computing the effects of refraction, would be expected to 
decrease the index of refraction and hence the calculated air mass. 

For comparison, values obtained by other investigators are also listed 
in Table 3. The values given by Treve appear to be the same,as those 
given in the Handbook of Geophysics and Space Environments, although 
no reference 1s given. The dif erence between the values calculated in 
this report and those of Treve ~ are most likely caused by the fact that 
Treve used the 1959 U.S. Standard Atmosphere, which differs slightly 
from the 1962 U.S. Standard Atmosphere above 20 km. 


The importance of the assumed profile fo, such refraction and air 
mass calculations has long been recognized.’ ~ Table 4 shows a compari- 
son of air mass for six different reference atmospheres. The 1962 U.S. 
Standard Atmosphere (58 lavers, 0 to 100 km) was taken from Reference 32, 
and the other five were taken from Reference 33. The air mass depends 
on wavelength, but this effect is small “ compared to the effect of 


different profiles. 
30 


Wilkes, M. V., "A Table of Chapman's Grazing Incidence Integral 
ch(X,x), Proc. Phys. Soc. London, B67, 304 (1954). 
31 


Valley, S. L. "Handbook of Geophysics and Space Environment," 

Pri 6: Standard Atmosphere 1962. 
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Edition) ,'"' Technical Report AFCRL-72-0497, Environmental Research 
Paper No. 411, 24 August 1972. 
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RESULTS 


It is a well known fact that the air mass increases rapidly with 
zenith angle for zenith angles greater than 90°. Fig. 3 shows quanti- 
tatively how rapidly the air mass increases for an observer at 45 km. 

It 1s seen that the air mass changes by 3 orders of magnitude over the 
range z = 90 to z = 96 degrees. Since the sun is not a point source, it 
can make a significant difference in the computed air mass whether a 
grazing ray comes from the top of the sun or from the bottom. This can 
be seen from Fig. 3 where the sun's finite angular diameter is shown 

for comparison. If the field of view of the spectrometer is constrained 
only to point somewhere on the sun, the resulting uncertainty in air 
mass can be large. This uncertainty is shown in Table 5 where the air 
mass and effective pressure are seen to vary by a factor of two over a 
range of 0.4° (the sun's angular diameter is about 0.5°). This effect 
has nothing to do with refraction, but is mentioned in passing because 
it can sometimes (depending on the experiment) be the major source of 
error, and errors caused by neglecting refractive effects must be com- 
pared with it. 


AS an example of the errors introduced by neglecting refraction, 
Fig. 4 was computed for the above example. The air mass was computed 
both including refraction (Fig. 3) and neglecting refraction 
(n = 1.0000000000). The ordinate of Fig. 4 is proportional to the 
percentage error, and the upper abscissa is the corresponding minimum 
height (including refraction). The observer altitude of 45 km was 
chosen as typical of the state of the art for high altitude balloon 
packages. It is seen in Fig. 4 that errors introduced by neglecting 
refraction (for the above geometry) are less than 10% when the minimum 
altitude is greater than 20 km. For minimum altitudes less than 20 km, 
the errors become rapidly larger. The result of neglecting refraction 
is to underestimate the minimum altitude, and therefore to overestimate 
the air mass. For a realistic atmosphere (where the density is a mono- 
tonic decreasing function with altitude), the minimum altitude computed 
from Equation 11 will always be less than the minimum altitude obtained 
by ray tracing. This causes the density at the minimum altitude (n_) to 
be overestimated, and this is the dominant term in the approximation to 
the Chapman function. 


Since Equation 11 shows the minimum height to be a function of 
observer height (as well as a number of other variables), one might 
expect that the error in minimum height would be a function of observer 
altitude. Fig. 5 shows this error for 3 different observer altitudes. 
Fig. 6 shows the corresponding errors in the computed air masses. It is 
thus seen that for a given grazing altitude, the above errors are larger 
for observers at higher altitudes. By extending this concept to an 
observer in a satellite at 200 km, the errors become quite significant 
as seen in Fig. 7. 


TABLE 5 


EXAMPLE OF POINTING ERROR 


OBSERVER ALTITUDE = 45 KM 


MIDLATITUDE WINTER MODEL 


ZENITH ANGLE AIR MASS 
95.8° 15.98 
96.0° 23049 
96.2” 30004 


EFFECTIVE PRESSURE (N/m*) 
14238.72 
20909.61 
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Figure 3. 
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Air mass vs zenith angle for an observer at 45 km. 
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Figure 4. Air mass error due to neglecting refraction for an 
observer at 45 km. 
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Figure 6. Error in air mass due to refraction for observers at 
20, 35, and 45 in. 
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Error in air mass due to refraction for an observer at 209 km. 
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Several investigators’?! have shown that the computed air mass 
depends on the atmospheric profile used. Fig. 8 shows the errors that 
would result for an observer at 45 km, if a winter profile were used to 
interpret summer measurements. Although the errors are not large (com- 
pared with, for example, pointing error), they can be easily reduced by 
using the atmospheric profile obtained on the day of the measurement. 


The sensitivity of the computed air mass with respect to wavelength 
errors was investigated by computing the air mass for an observer at 45 
km, at both 5 and 0.5 micrometers. Even for the minimum altitude of 
less than 10 km, the difference in computed air mass was less than 0.4%, 
and the difference in minimum altitude was 19.4 meters. That the air 
mass is not a sensitive function to small changes in refractive index is 
consistent with the neglect of the contribution of water vapor and 
anomalous dispersion to the index of refraction. 


Differential petiwetimac*"* which has been studied by several in- 
vestigators, has not been included here. It is expected to be small for 
balloon altitudes, but might possibly be significant at satellite alti- 
tudes, 


SUMMARY 


Effects due to refraction have been investigated for several repre- 
sentative occultation geometries. Errors in the computed air mass and 
effective pressure, associated with neglecting refraction are shown to 
be a function of observer altitude. For observers in high-altitude 
aircraft, the errors introduced by neglecting refraction are probably 
less than those introduced by uncertainties in solar pointing. For 
observations from high altitude balloons or satellites, refraction 
should definitely be taken into account when the minimum altitude is 
less than 20 km. For all cases, accurate knowledge of the astronomical 
zenith angle is a necessity. 


For interpretation of infrared transmission spectra, where the 
absorption 1s pressure dependent, errors in minimum altitude are doubly 
important. Not only will the air mass be overestimated but also the 
effective pressure will be overestimated. By neglecting refraction, the 
mixing ratios for minor constituents in the lower stratosphere would 
tend to be underestimated. 
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APPENDIX A. AIR MASS TABLES 


In order to present the results of this research in a form that 
would be convenient for other investigators, a series of refraction air 
mass tables was generated. Two atmospheric models were assumed: a 
midlatitude summer and a midlatitude winter. Other models could have 
easily been included, but these two were judged adequate for most appli- 
cations. 


The temperatyze profiles and the surface pressures are those of 
McClatchey et al. The pressure profile is calculated from the temper- 
ature profile and the surface pressure by integrating the hydrostatic 
equation. These calculated pressures agree to within 2% of those of 
McClatchey, up to 30 km. The layer thickness is 250 m up to 40 km, 

1000 m up to 70 km, and 5000 m up to 100 km. There are thus 197 layers 
in these model atmospheres. 


Tables are given for observer altitudes between 10 and 50 km, in 
steps of 1 km. For each observer altitude, the zenith angles range from 
80° to 89° in steps of 1°, and from 89° to the solid earth in steps of 
0.2°. A single wavelength of 5.0 micrometers was chosen, since the air 
mass is not strongly wavelength dependent. 


The first column, labeled APP. Z, represents the angle of arrival 
(deg) with respect to the local normal. Columns 2 and 3 represent the 
observer height (m) and the highest level (m) included in the atmos- 
pheric profile respectively. The minimum height (m) in column 4 repre- 
sents the observer altitude for z less than or equal to 90°, and the 
tangent altitude for a grazing ray, for z greater than 90°. Columns 5 
and 6 give the air mass from the Smith and Smith approximation to the 
Chapman function, and the ray tracing numerical integration, respec- 
tiygly. In thesg two columns, 1.0 air masses is defined as 2.15335 x 
10°" molecules/m'. HDIFF (m) represents the difference between the ray 
traced minimum altitude, and that obtained from Equation (11), neglecting 
refraction. Columns 8 and 9 give the refraction (deg) and the astronomical 
zenith angle (deg). Columns 1 gnd 8 add to give column 9. The single 
layer effective pressure P (N/m’) and temperature T (°K) for a uniformly 
mixed gas are given in columns 10 and 11 respectively. Golumn 12 gives 
the pressure at the tangent or minimum altitude P ean (N/m As 
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APPENDIX B. SOURCE STATEMENTS FOR THE RAY TRACING PROGRAM 


The following listing is for the FORTRAN IV program used to cal- 
culate the values in Table 5. The calculations were performed on the 
BRLESC II computer, at the Ballistic Research Laboratories, which has 
a Single precision accuracy of 16 significant digits. For use on other 
computers which carry fewer significant digits, the use of double pre- 
cisions should be investigated. 


The data deck consists of the following cards: 
Card 1. Title (13A6) 
Card 2. LEVELS (15): the number of levels in the 


atmospheric profile to be read in. 


ALTORH (F5.0): variable defining input altitudes 
as geometric (#0) or geopotential (=0) meters. 


ALAT (F10.8): latitude of the observer in degrees. 
NOPT (110): format variable for the input profile. 
Card 3 to LEVELS + 2 input profile, altitude, temperature, 
and pressure (units depend on NOPT) 
Card LEVELS + 3: data card for specific geometry 
DD1 (F10.7) = zenith angle (degrees) 
DD2 (F9.6) = observer height (meters) 


DD3 (F10.6) = source height = highest level of 
profile (meters) 


DD4 (F11.3) = wavelength (micrometers ) 
STOP (14): program option parameter 


0 = Apparent zenith angle given - 
compute astronomical 

1 = Astronomical zenith angle given - 
compute apparent 

2 = calculate preceding case with 
new wavelength 

3 = call exit 


ALTORH (F8.1): parameter for geometric (#0) or 
geopotential (=0) meters for observer height. 


Any number of cards for specific geometries can be used (within time 
limit). To end the computation, the last data card should have a ''3" 
in column 44. 
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